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The reformulated coupled-cluster method (CCM), in which average many-body potentials are 
introduced, provides a useful framework to organize numerous terms appearing in CCM equations, 
which enables us to clarify the structure of the CCM theory and physical importance of various 
terms more easily. We explicitly apply this framework to 4 He, retaining one-body and two-body 
correlations as the first illustrating attempt. Numerical results with using two modern nucleon- 
nucleon interactions (AV18 and CD-Bonn) and their low-momentum interactions are presented. 
The characters of short-range and many-body correlations are discussed. Although not considered 
explicitly, the expression of the ground-state energy in the presence of a three-nucleon force is given. 

PACS numbers: 21.10.Dr, 21.60.-n, 31.15.bw 
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I. INTRODUCTION 



The understanding of various properties of atomic nu- 
clei starting from the realistic nucleon-nucleon interac- 
tions is one of fundamental problems in the theoretical 
nuclear physics. The recent studies in this field are char- 
acterized by several new developments. Resulting from 
developments through 1960's and 1990's, there are now 
various parameterizations of the nucleon-nucleon inter- 
action which reproduce experimental two-body scatter- 
ing and deuteron data with high precision. Systematic 
introduction of a three-nucleon force is in progress espe- 
cially in the framework of the new development of the 
potential description in the chiral effective field theory 
[ll, 0- Progress is also seen in an effective interaction 
theory in a restricted space, driven by active studies of 
low-momentum interactions Q. In addition, various ab- 
initio frameworks of many-body calculations have been 
explicitly applied to atomic nuclei, such as a Monte-Carlo 
method [|J, a no-core shell model Q and a coupled- 
cluster method (CCM) @. Among them, the CCM is 
a promising method toward heavier nuclei because of the 
advantage of the size-extensivity. 

The CCM was devised for many-body problems in nu- 
clear physics in 1950's @, H| and the achievements in 
1970's by the Bochum group were reported in Ref. @. It 
was, however, almost discarded in nuclear physics com- 
munity. The CCM found its place in quantum chemistry 
[ToL [Tl| as a tool of the first-principle calculation. The 
method was reintroduced around 2000 in the description 
of atomic nuclei @, UH . Because there are specific dif- 
ficulties in a description and a treatment of the short- 
range part of a nucleon-nucleon interaction, applications 
of the CCM to nuclei require deeper understanding of an 



effective interaction theory. It is not surprising that a 
renewed interest arises in the CCM in parallel with the 
development of low- momentum interaction theory [3[ . 

There is already a considerable number of CCM cal- 
culations of light nuclei in recent years [1 S 17 1. These 
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include the extension to excited states 
far from the stability line [15| and heavier nuclei such as 
40,48 Ca [3, [l||. Contributions of a three-nucleon force 
has also been estimated in this method [l8j]. Neverthe- 
less, because of the usefulness of the CCM framework 
to solve the quantum many-body problem as exactly as 
possible, it is important to clarify the structure of this 
framework as much as possible and to try to obtain the 
physical understanding of various correlations in a trans- 
parent way. Such an attempt was proposed by Suzuki 
[19| in the early 1990's as an application of the similarity- 
transformation theory for a quantum-mechanical eigen- 
value problem. His reformulation introduces average 
many-body potentials as the generalization of the one- 
body mean field. The concept of the mean field is essen- 
tially important in almost all many-body systems, as the 
independent single-particle model is empirically estab- 
lished in those systems. As is shown below, the many- 
body average field is useful to group various terms which 
appear in the CCM practical calculations, and thus is 
helpful to elucidate the structure of the CCM many-body 
correlations. The transformation of the Hamiltonian to 
the normal ordered form with respect to the reference 
state is inherent in this formulation. Unfortunately, ex- 
plicit applications of this method to nuclei has not been 
carried out so far. It is useful at the present stage to pro- 
mote this formulation for the description of actual nuclei 
with using the realistic nucleon-nucleon interaction, in 
view of the future wide use of the CCM framework in 
nuclear many-body problems. 

We recapitulate the CCM formulation with introduc- 
ing average many-body potentials in Sec. II, starting 
with the basic idea of the CCM. Although we do not 
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consider the contribution of a three-nucleon force in the 
present numerical calculations, the expression including 
three-body terms are presented in the initial part of Sec. 
II. As the first attempt of the actual application of this 
method to nuclei, we employ the approximation in which 
only one-body and two-body amplitudes are retained. 
This truncation is referred to as CCSD in the literature. 
Numerical calculations are carried out for 4 He with us- 
ing two modern nucleon-nucleon potentials, the CD-Bonn 
potential [2C§ and the Argonne AV18 potential [21|. Re- 
sults of the energy and radius of 4 He in a harmonic oscil- 
lator model space are presented in Sec. III. First, results 
with the bare potentials are shown, and then the calcu- 
lations with low-momentum equivalent interactions are 
discussed. The magnitude of the contributions from the 
average many-body potential is estimated. The indepen- 
dence of the results on the harmonic oscillator constant of 
the reference state is demonstrated. If we do not include 
one-body amplitudes, the results varies with changing 
the oscillator constant. Conclusions follow in Sec. IV. 

II. COUPLED-CLUSTER METHOD 

We first give the idea of the CCM in its naive form, 
and next present the introduction of many-body average 
potentials. 

A. Basic 

The basic ansats of the coupled-cluster method @, HI 
is that the exact state |^) of the Hamiltonian H of A 
particles with the mass m is given by the transformation 
operator e s acting on the reference state |$o) : 

l*)=e s |<fo). (1) 
The Hamiltonian H consists of the one-body kinetic en- 
ergy ti — 5^rfejj the two-nucleon interaction Vij, the 
three-nucleon interaction Uy^, and so on. To be explicit, 
we write the Hamiltonian as follows. 

H = E ^ -Tc+Y, «« + E ( 2 ) 

i i<j i<j<k 

where Tq denotes the center-of-mass kinetic operator. 
In the following, we absorb the one-body part of Tq in 
the one-body kinetic energy, namely the operator being 
ti = ^^-■^k 2 , and the two-body part of Tq in the two- 
body interaction . 

Solving the time-independent Schrodinger equation 

He s \^ ) =£ e S |$o), (3) 

is to determine the transformation operator S, which con- 
sists of one-body part, two-body part, and up to A-body 
part. 

s = y] Si + y] + s i:jk -\ — hsi2...A- (4) 

i i<j i<j<k 



We refer to single-particle states in |$ ) as hole 
states denoted by h and other single-particle states 
as particle states denoted by p. By definition, 
only the particle (hole) states appear in the bra 
(ket) configuration of the transformation amplitude 
(piP2- ■ ■Pn\si2--- n \hih 2 - ■ ■ h n )A: where the suffix A 
means an antisymmetrized matrix element; namely, 
\h 1 h 2 )A = \hih 2 - h 2 hx), \hh'h") A = \hh'h" - hh"h' + 
h'h"h - h'hh" + ti'hti ~ h"h'h) and so on. 

The prescription of the CCM is to rewrite the 
Schrodinger equation © as 

e- s He s \<$> ) = £ |$o>- (5) 

It is noted that in the standard CCM theory, the Hamil- 
tonian is supposed to be normal-ordered with respect to 
the reference state at this stage. Hence we denote the 
Hamiltonian by H n . The operator S is determined by 
the decoupling equation for the similarity-transformed 
Hamiltonian e~ H n e s : 

{np-nh\e~ s H n e s \^) = 0, (6) 

with n — 1,2, ■■ -A. Here, (np - nh\ represents an ar- 
bitrary rt-particle and n-holc state. Determining 5, the 
energy is given by 

E Q = <$ |e- S -ff„e s |$o). (7) 

The similarity-transformation e~ s H n e s may be eval- 
uated systematically by the Baker-Campbell-Hausdorff 
formula: 

e- s H n e s = H n + [H n , S] + i [[H n , S], S] + ■ ■ ■ . (8) 

In practice, the Hamiltonian should be rearranged in the 
normal-order form with respect to | <I>o) The trans- 

parent procedure of the transformation in view of the ef- 
fective interaction theory was formulated by Suzuki [l9| . 
This method with introducing many-body average po- 
tentials is outlined in the next subsection. 



B. Many-body average potential 

When the original Hamiltonian H do not have 
more than three-nucleon interactions, the expansion of 
e~ s He s terminates with the four-folded commutator and 
produces various terms up to six-body operators. If 
a three-nucleon interaction is present, more than five- 
folded commutators appear. Suzuki and collaborators 
(l9l |23| reformulated the CCM by introducing many- 
body average fields, which enables us to write down the 
decoupling equations in a compact form and thus to clar- 
ify the structure of various many-body correlations. We 
manipulate explicit expressions in this formulation and 
carry out numerical calculations. 

The similarity-transformed Hamiltonian may be reor- 
ganized by introducing an auxiliary many-body potential 

U = E Ui + E Ui i + E Ui i k + " ■ • ( 9 ) 

i i<j i<j<k 
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FIG. 1: Bubble insertion diagram for Hartree-Fock mean field. 
The antisymmetrization is not shown explicitly. 



i i<j 

+ X (^' fe _ ^ — 

i<j<k 



(11) 



The prototype of the auxiliary potential for a many-body 
system is a Hartree-Fock one. The one-body Hartree- 
Fock mean field for states \a) and |6) is defined by folding 
two-body interactions with respect to occupied states \h): 



The transformed Hamiltonian is written as 
e - s He s = e~ s (H + U)e s - e~ s Ue s 

i i<j i<j<k 



+ •••} 



■}.(10) 



i<j i<j<k 



That is, 



e S {J2 t i + J2 v v + X v ^ eS 

i i<j i<j<k 



\u H F\b) = ^(ah\v 12 \bh) J 



(12) 



The summation is depicted as a bubble insertion diag ram 
of Fig. 1. The formulation by Suzuki et al. |l9l [23| 
generalize the HF potential to many-body states. That 
is, the auxiliary potential U is determined so as to cancel 
all the bubble-insertion contributions. 



(-l) fe+1 

ai ■ ■ ■ a n \ui... n \bi ■ ■ ■ b n ) A = 2J rj X ( fll " ' anhl ' ' ' hk\vi... n +k\h ■ ■ ■ b n h\ ■ ■ ■ h k ) Al 



k>l 



(13) 
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where a and b stand for either particle or hole state. It is instructive to give explicit expressions for the one-body u\ 
and the two-body w 12 . 



ti| 6) = ^2(ah\vi 2 \bh) A - ^^(ahh'\vi 2 3\bhh'} A + ^^2(aia 2 hh'\vi23i\bib 2 hh') 



A + 



h W 



hh' 



iia 2 \u 12 \bib 2 ) a = ^2(aia 2 h\vi 23 \bib 2 h) a - Try y^ y (aia 2 hh'\vi 234: \bib 2 hh') 



(14) 



(15) 



hh' 



It is apparent that u is a sum of bubble-insertion contributions of the similarity-transformed interactions, if we rewrite 
the above definition of u\ as 



(a|t*i|6) = '^(ah\v 12 - u 12 \bh) A + ^(ah\u 12 \bh) A - -^^(ahh'\v 123 \bhti 



(16) 



- J2 {ahh'h"\v l234 \bhh'h") 



hh' 



3! 



hh'h" 



'^{ah\vi 2 - ui 2 \bh) A + ^ y^^ahh! \vx23\bhh!) a + (-L - ^ J X { ahh ' h"\vi 23i \bhh' h") 

h hh' \ ' ■ / hhh „ 

X( aft -I«i2 - ui 2 \bh) A + i ^2(ahh'\vi23 - ui 23 \bhh') A + ^ X {ahh' h" \vx234\bhh' h") a 



fc.fr/ hhh" 



(17) 



and also for mi 2 

(aia 2 |-u 12 |6ifo 2 )A = X^ aia2/l l^ 123 - u w \bxb 2 h) A + (a^hlu^b^h) A } - ^ (ai a 2 fefo' 1 5i234 1 &i &2 hh') a H 



/!/,' 
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= ^(a 1 a 2 h\v 1 23-u 1 23\b 1 b 2 h)A + ^^(a^hh'lv^^b^hh^A^ ■ (18) 



2 

hh 



Therefore, e~ s (H + U)e s is now arranged in the normal- 
ordered form with respect to the reference state. It means 
that the following decoupling conditions should be im- 
posed. 

(np-nh\e' s (H + U)e s \<S> ) = 0. (19) 

This corresponds to considering the decoupling condition 
for the normal-ordered Hamiltonian H n in the standard 
CCM, Eq. ©. The potential U thus defined may be 
called a many-body average potential. 

To proceed to actual calculations, explicit expressions 
of the similarity-transformed interaction e~ s He s are 
needed. The result of the straightforward calculation is 
given in the Appendix. In these expressions the terms 
including a three-nucleon interaction and a three-body 
amplitude Sijk are retained. However, because we do not 
consider these contributions in the CCSD truncation in 
this paper, we do not keep these terms here after. 



interactions V12 and the correlation amplitudes s± and 
si 2 . Because hi — e~ Sl tie Sl +iti = (1 — si)ii(l + si) + ui, 
the one-particle-one-hole decoupling equation becomes 



{pfh^h} = (p\(l-s 1 )t 1 (l+s 1 ) + u 1 \h)=0. (20) 



The explicit expression of (p\u\\h) is given below, Eq. 
(|2~7| . Next, the two-particle- two- hole decoupling equa- 
tion {pip 2 \vi_ 2 \h\h 2 ) a — reads 

= (pip 2 |(l -31-82 + 8iS 2 - Si 2 )fl2 

x(l + si + s 2 + sis 2 + si 2 ) + Ui2\hih-2)A = 0. (21) 



C. Decoupling equations and energy 

Figuring out the explicit form of the matrix element 
Some manipulation is necessary to obtain explicit ex- (pip 2 \ui 2 \h\h2) a according to Eq. (|18|) . this equation 
pressions of matrix elements of u in terms of the original can be written in the following rather compact manner: 



where 



(P1P2KI - si - s 2 + sis 2 + s 12 )ui2(l + si + s 2 + s 1 s 2 + si 2 )\hih 2 )A 
+(j>ip 2 \(ei + e 2 )si2 - si 2 (ei + e 2 )|fti/i 2 )A + {piP2\uf 2 + uf 2 \h 1 h 2 }A = 0, 



(hi\ei\h 2 ) = (hilhil + s^+U!^), 
(pi\ei\p 2 ) = (pi|(l-si)ti+ui|p 2 ), 

(hi\ui\h 2 ) = y~]{hih'\vi2(l + si + s 2 + S1S2 + si 2 )\h 2 h')A, 



(Pl|ui|p 2 ) = y]{pih'\(l - Sl)«l 2 (l + Si + S 2 )\p2h')A - ^(PIP'ISUVU&P') 



(pip 2 \uf 2 \hih 2 ) a 



A- 



(22) 



(23) 
(24) 

(25) 
(26) 



(p\ui\h) = (ph'\t 2 si2 + (1 - si)wi 2 (l + si + S2 + sis 2 + s 12 )\hh') A 

+ \ Y,i 2 ( h ' h "\ V ^P' h ")A{PP'\-S 1 2\hh') A + {h'h"\v l 2\p'h) A {pp'\s l 2\h'h")A} 

h'h" p' 

+ \ E \ E 2(h'h' / \v 12 \py)A{2(pp'\ Sl 2\hh')A(p , '\s 1 \h") + (pp'\si2\h'h") A (p"\si\h)}, (27) 

h'h" p'p" 



E 

h'p' 



{(P2h'\(l - si)vi2\p'hi) A + ^T(M'|(1 - si)vi 2 \p'p") A (p"\si\hi)}{p lP '\s 12 \h 2 h') A 



-{{p 2 h'\{l - si)vi 2 \p' h 2 ) A + y^(p 2 ^l(l - si)vi 2 \p'p") A {p"\si\h2)}(pip'\si2\hiti) A 

p" 



■5 



-{(piti\(l - sxMp'h^A + X>i^'l(l - s 1 )v 12 \ P 'p") A (p"\ Sl \h 1 )}{ P ' P2 \s 12 \h 2 h') 



p" 



-{(pih'\(l - si)vi 2 \p'h 2 ) A + ^(pih'\(l - s 1 )v 12 \p'p") A (p"\s 1 \h 2 )}(p'p 2 \s 12 \h 1 ti) 



( Pl p 2 \uf 2 \h 1 h 2 )A = ±\ ^(^^>12bV%{Wk 12 |/l 1 / l %(p 2P '>12|/i 2 /l" 
h'h" p'p" 



-{P\p'\si2\h 2 h') A (p 2 p"\s 12 \h 1 h") A }. 



, (28) 



(29) 



Vl2 



Vl2 




N\ M 

M \fa P2 / M 



FIG. 2: Diagram representation of the S12 contribution to 
the one-body average potential ui, Eqs. (|25|) and (|26|l The 
antisymmetrization is not shown explicitly. 




hi h 2 pi p2 




FIG. 3: Diagram representation of the S12 contribution to the 
two-body average potentials uf- 2 and Ui 2 , Eqs. (|28p and (|29p . 
respectively. The antisymmetrization is not shown explicitly. 



D. Energy and root- mean-square radius 

The ground-state energy is obtained by Eq. (J7]). In the 
case that only two-body interactions are present, ampli- 
tudes higher than three-body Sijk do not contribute to 
the CCM energy and the expression becomes: 

E Q = J2{h\ti(l + Sl )\h) + ^J2( hh '\ v ^ 

h hh' 

x(l + si + s 2 + s x s 2 + s 12 )\hh') A . (30) 

Effects of the higher order correlations enter through the 
one-body and two-body decoupling conditions. Although 
we do not consider contributions of a three-nucleon force 
in this article, it is instructive to present the energy ex- 
pression in the presence of and Syfc, which is given in 
Appendix D. There we further rewrite the expression by 
introducing an effective two-nucleon interaction obtained 
from a three-nucleon interaction. 

In the CCM framework, the evaluation of expecta- 
tion values (vflQI^) of other observables Q than the en- 
ergy needs explicit expansion of e s Qe s . Although only 
linked-cluster diagrams are needed to be taken into ac- 
count, a perturbative estimation in all-orders is not fea- 
sible. If the energy is determined variationally, we may 
use the Feynmann-Hellman theorem to obtain the matrix 
element of the given observable Q by solving the problem 
of the constrained Hamiltonian H + XQ. That is, sup- 
posing the state |^(A)) to be the solution of the forced 
Hamiltonian, the following relation holds. 



Note that £ |W*'MW»'I (i E PP ' W)a(pp'\) is a 
unit operator in an occupied (unoccupied) space. The 
S12 contribution to the single-particle potential ui, Eqs. 
(|2"5|) and (|26l) , is diagrammatically represented as in Fig. 
2, and the potentials uf 2 and uf 2 , Eqs. (I2"5t and (|2U|) . in 
Fig. 3. It can be checked that Eq. (|2ip completely agrees 
with the CCSD decoupling equation in the literature, as 
it should: e.g., Eq. (35) of ref. @. 

We need more explicit expressions to evaluate matrix 
elements such as (pip 2 \(si + s 2 )vi 2 \hih 2 ) A , and algebraic 
calculations of the angular-momentum coupling for the 
numerical calculations of these expressions. We solve 
these non-linear equations for Sj and Sij in an iterative 
way, using essentially the Newton method. 



(*|Q|*)= hm— (*(\)\H 



AQ|*(A)>. 



(31) 



Because the CCSD energy does not have a rigorous vari- 
ational property, (^ , (A)|i?|^(A)) may depend linearly on 
A, which introduces some uncertainty to the calculated 
expectation value of Q. Nevertheless, we calculate the 
root-mean-square (r.m.s.) radius of the 4 He by this 
method, which is a square root of the expectation value 
of the operator 



Qr 



1 A 
IB* 



Rf 



.4 



A 2 



(32) 
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where R = 4 J^i r i is the center-of-mass coordinate. Be- 
cause the CCSD energy reaches almost variational limit 
in practice, the estimation by the constrained calcula- 
tion is more reliable than the perturbative estimation in 
low-orders. As is shown explicitly in the next section, 
the reliability of the prescription of Eq. (j3Tj) is suggested 
by the fact that the calculated r.m.s. radii are almost 
independent on the oscillator constant of the harmonic 
oscillator basis. 

It is helpful to note that the r.m.s. radius of the 4 He 
reference state |$o) = |(0s) 4 ) with the oscillator constant 
v is given by 

vWi) = -^j=- (33) 



III. RESULTS 

We carry out numerical calculations in the harmonic 
oscillator basis. The oscillator constant is denoted by v 
which is related to the frequency ui as hoj = — v with 
m being the nucleon mass. The size of the harmonic- 
oscillator basis is specified by the quantum number 
N m ax — {2n+i) max , where n is a nodal quantum number 
[n — 0, 1, 2 • ■ ■) and I an orbital angular momentum. The 
number of the major shells included is given by N max + 1. 
To see the dependence of CCM calculations on the size 
of the model space, we use four different choices of N max , 
as shown in Fig. 4. 

As the first application of the reformulated CCM, we 
consider only 4 He, for which the reference state |$) is 
taken to be (0s) 4 . To keep the computaional cost low, 
we set N max < 22 and introduce the cutoff at the upper 
right corner of the model space as in Fig. 4. Two-body 
partial waves up to J — 6 are included. We employ 
two realistic nucleon- nucleon interactions; the CD-Bonn 
potential 0| and the AV18 potential [2l|. The charge 
dependence of the CD-Bonn potential is averaged in ac- 
tual calculations for the use of the isospin basis. These 
two potentials arc different in their strength of the ten- 
sor force. The former is known to have relatively weak 
tensor components. Discussions are chiefly presented for 
the results of the CD-Bonn potential. 

Two-body matrix elements of the bare nucleon- nucleon 
interaction in the harmonic oscillator basis are evalu- 
ated by the Talmi-Moshinsky transformation. Because 
the computation of Talmi-Moshinsky coefficients is time 
consuming when the oscillator quantum number becomes 
large, we use the approximation that the center-of mass 
(CM) quantum number is restricted to be (2n+£)cM < 6. 
We have checked that the error due to this simplification 
is negligibly small. 




10 



6 



6 10 14 18 22 
2n 1+ li 

FIG. 4: Harmonic-oscillator model space for CCM. Four 
choices, (2n + £) max — 6, 10, 14, and 22, are used for numerical 
calculations. 



A. Calculations with CD-Bonn bare interaction 

Calculated CCSD ground-state energies of 4 He with 
the CD-Bonn potential are shown in Fig. 5, as a function 
of the size of the harmonic-oscillator base, N = 2n + £, 
with the oscillator constant v — 0.56 fm~ 2 . As a ref- 
erence, the energy of the Fadeev-Yakbovsky calculation 
(2~4| with the same interaction is also shown. Increas- 
ing the size of the model space, the CCSD energy goes 
down toward to the F-Y energy. However, at N rnax = 22, 
the energy is still short by about 8 MeV to the reference 
value of the Faddeev-Yakubovsky calculation. When the 
AV18 potential is employed, the shortage is about 21 
MeV, compared with the corresponding F-Y result of - 
24.3 MeV. This result is not surprising, because the s.p. 
space is not sufficient to properly treat the short-range 
correlation of the ordinary realistic nucleon-nucleon in- 
teraction. The momentum scale in the model space may 
be estimated by the expectation value of the square of 
the momentum operator (p 2 ) = (N, nax + |)za The num- 
ber is yj (p 2 ) — 3.6 fm _1 for N max — 22 and v = 0.56 
fm~ 2 , which is much smaller than the high momentum 
relevant to the nucleon-nucleon short-range correlation. 

It is difficult in practice to employ a large harmonic 
oscillator basis for CCM calculations enough to deal 
with the nucleon-nucleon short-range correlation. The 
established wisdom in nuclear physics is to treat the 
short range singularity by solving the G-matrix equa- 
tion as a two-body problem in the nuclear medium 
[25L [26| . Though the perturbative expansion in terms 
of G-matrices is a solid framework, the CCM provides a 
more systematic and compact way to treat higher-order 
correlations. One standard approximation besides the G- 
matrix approach is to use a two-step method, in which 
two-body correlations are initially treated in a large space 
and then solve many-body correlations in a restricted 
space. This strategy is adopted, e.g., in the unitary- 
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FIG. 5: Dependence of CCSD ground-state energies of 4 He 
on the harmonic-oscillator model-space size N ma x with the 
bare CD-Bonn potential [2(j. The oscillator constant is v — 
0.56 fm" 2 (huj = 23.2 MeV). The energy by the Faddeev- 
Yakubovsky calculation [24| is shown for reference. 



FIG. 6: Dependence of CCSD ground-state energies of 4 He on 
the harmonic-oscillator model-space size (2n+£) m ax with low- 
momentum equivalent interactions of three cutoff momenta, 
A = 4.0, 3.0, and 1.9 fm -1 . The bare interaction is the CD- 
Bonn potential and the oscillator constant is v — 0.56 fm -2 
(huj = 23.2 MeV). 



model-operator approach [27j |. 

A different method has been recently developed. 
Namely, the bare interaction is transformed to the 
half-on-shell two-body equivalent interaction in low- 
momentum space 0, which is free from the high- 
momentum singularity. This interaction has been 
shown to be used in perturbation calculations in a low- 
momentum space. Though the linkage between the har- 
monic oscillator basis and the low-momentum space is 
somewhat obscure and therefore the treatment in the har- 
monic oscillator basis from the beginning is desirable, we 
use the low-momentum space interaction in this paper 
because of its simplicity and possibilities of the compar- 
ison with calculations in the literature flij . 

It is noted that in a restricted model space more than 
three-body correlations involving high-momentum states 
are totally dropped. These contributions are to be recov- 
ered by more than three-nucleon induced interactions in 
a restricted space. Otherwise some adjustable param- 
eters may be introduced, if we take a viewpoint that 
high-momentum components of bare nucleon-nucleon in- 
teractions are themselves cannot be determined without 
uncertainties. 



B. Calculations with low-momentum interaction 

In this subsection, calculated results for low- 
momentum space equivalent interactions of three cutoff 
momenta, A = 4.0, 3.0, and 1.9 fm -1 , are presented. The 
CCSD 4 He ground-state energies and radii with the CD- 
Bonn potential are shown in Figs. 6 and 7 as a function 
of the size of the harmonic-oscillator basis. 

Results with A = 1.9 fm^ 1 show that the energy of 
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FIG. 7: Same as in Fig. 6, but for CCSD ground-state point- 
nucleon r.m.s. radii of 4 He. The experimental point proton 
radius is shown as a guide, which is taken from ref. |28f ] . 



the Faddeef-Yakubovsky calculation with the same low- 
momentum equivalent interaction is almost reproduced 
already at N max = 6. This implies that if the two-body 
correlations are renormalized at A = 1.9 fm^ 1 , more 
than three-body correlations in the model space scarcely 
change energy. The evaluated radius is also insensitive 
to the size of the model space for A = 1.9 fm -1 . On 
the other hand, the difference between the energy with 
A = 1.9 fm -1 and the F-Y energy with the bare inter- 
action suggests that more than three-body correlations 
involving higher- momentum components are of the order 
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FIG. 8: Oscillator constant dependence of the ground-state 
energies of 4 He for the model space N ma x = 10 with using the 
low- momentum equivalent interaction of A = 3 fm _1 . The 
solid curve is the result of full CCSD calculation. The dashed 
curve shows the result with discarding uf^ and itf 2 . The dot- 
ted and dot-dashed curves are results of the calculation for 
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with and without w^ 2 and wf 2 , respectively. 



FIG. 10: Oscillator constant dependence of the ground-state 
energies of 4 He for the model space N max — 22 with using the 
low- momentum equivalent interaction of A = 3 fm _1 . The 
solid curve is the result of full CCSD calculation. The dashed 
curve shows the result with discarding u^ 2 and uf 2 . The dot- 
ted and dot-dashed curves are results of the calculation for 

,B 
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with and without u^ 2 and wf 2 , respectively. 



n 1.61- 

CO 



^ 1.4h 

CD 



1.2- 
15 




20 25 30 

ha>[MeV] 



35 



n 1.6 

CO 



- 1.4 

CD 



1.2 



\\\ 



AV18 



3T- 





\ W CD-Bonn 



15 20 25 30 

hco[MeV] 



35 



FIG. 9: Same as in Fig. 8, but for ground-state point-nucleon 
r.m.s. radii of 4 He. The thin dot-dashed curve is almost 
indistinguishable from the bold dot-dashed curve. The r.m.s. 
radius of the reference state |$o), Eq. (|33[) , i s shown by a 
two-dot-dashed curve. 



FIG. 11: Same as in Fig. 10, but for ground-state point- 
nucleon r.m.s. radii of 4 He. The thin dot-dashed curve is 
almost indistinguishable from the bold dot-dashed curve. The 
r.m.s. radius of the reference state |$o)> Eq. (|33|) . is shown 
by a two-dot-dashed curve. 



of a few MeV. Observing that the energy at A = 3 fm _1 
coincides with the F-Y energy, we see that the three-body 
correlations at the momentum range of 2 ~ 3 MeV/c 
is not negligible. It is remarked that the energy with 
A = 1.9 fm _1 is slightly lower than the F-Y energy with 
the corresponding low-momentum interaction. This fea- 



ture is different from the tendency in ref. 

Next, we turn to the oscillator constant dependence 
of the calculated results for CD-Bonn and AV18. This 



case, we also show results of the calculation in which S\ 
is set to be and/or the potentials uf^ in Eq. (|2"Tj) arc 
omitted. Figs. 8 and 9 are the results for the model space 



N„ 



10 , and Figs. 10 and 11 for N„ 



22 with 



using the low-momentum equivalent interaction of A = 3 
fm _1 . Results from the CD-Bonn (AV18) potential are 
shown by thick (thin) curves. 

The full CCSD results are seen to be almost inde- 
pendent on the oscillator constant both for ground-state 
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energies and radii of 4 He, except for the region below 
fiu ~ 20 for N max — 10. The tendency in the case 
of N max — 10 that the energy with larger fku pro- 
vides slightly lower energy is understandable, because the 
s.p. wave functions of the larger oscillator constant con- 
tain relatively larger high-momentum components. This 
point, however, needs to be clarified in detail together 
with the slight convex shape of the oscillator constant 
dependence of the energy. The difference of the solid 
and the long-dashed curves in Figs. 8~11 shows that the 

A B 

neglect of the potentials u 12 ' does not much affect the 
energy and the radius. The ground-state energy becomes 
less attractive by about 0.5 MeV and the radius becomes 
slightly smaller by about 0.02 fm in the both cases of 
CD-Bonn and AV18. 

The difference of the contributions of correlations de- 

A B 

scribed by u 12 ' with CD-Bonn and AV18 reflects the dif- 
ferent character of these interactions, typically observed 
in the strength of the tensor components. Together with 
the convergence with respect to the size of the model 
space, we see that the many-body problem with the low- 
momentum equivalent interaction is solved adequately by 
the CCM framework at the CCSD level. 

If we neglect s\ terms from the beggining, the CCD 
approximation, calculated energies and radii vary with 
changing the oscillator constant of the model space, as 
Figs. 8~11 show. Increasing the frequency huj, the en- 
ergy attains a minimum value at hui = 21 ~ 23 MeV 
and the radius decreases monotonically in almost paral- 
lel with the r.m.s. radius of the reference state \$o), Eq. 

Even in this case, the effect of it 12 is comparable 
to the case of fully including s%. It is interesting to ob- 
serve that the difference between the energy with s% = 
at its minimum value and the full CCSD energy is rather 
small. The r.m.s. radii with and without s\ are also seen 
to coincide each other at the point where energies with 
and without Si become close. 

The order of 0.5 MeV of the difference between calcu- 
lated energies with and without the uf^ B terms is smaller 
than that of the energy gain due to the two-body high- 
momentum correlations. This implies that the two-body 
correlations including high-momentum components are 
primarily important and the remaining many-body cor- 
relations summed up by uf^ B in a restricted model space 
plays a minor role in the magnitude of energy. This is 
consistent with the conventional understanding, regard- 
ing the justification of the independent particle picture 
in spite of the short-range repulsion of the bare nucleon- 
nucleon interaction. 

The above observation suggests that three-body corre- 
lations including high-momentum components are worth 
to investigate quantitatively. On the other hand, In the 
renormalization-group view-point on which the develop- 
ment of low-momentum interaction Q is conceptually 
based, such correlations should be hidden in adjustable 
parameters as far as low-energy properties of nuclei are 
concerned. When the experimental data at intermediate 
energy in nuclear physics is accumulated in the future, 



these correlations will be addressed as the important sub- 
ject to analyze. 



IV. SUMMARY 

We have evaluated the ground-state energy and ra- 
dius of 4 He in the reformulated coupled-cluster method 
in which many-body average potentials are introduced. 
This method was proposed by Suzuki in 1990's [H, HII, 
but has not been practiced for actual nuclei. The formu- 
lation clarifies the procedure to organize the Hamiltonian 
in a normal-ordered form with respect to the reference 
state, and thus the structure of the CCM approach as 
the many-body theory. The obtained decoupling equa- 
tion in the truncation at the two-body correlation level, 
commonly referred to as the CCSD approximation, is ar- 
ranged in a concise form as given in Eq. (|22|) by virtue of 
the average potential, although the content is necessarily 
not different from that of other CCM calculations Q. 

Numerical calculations are carried out in the harmonic 
oscillator basis. Major shells up to N max = 22 are pre- 
pared, while the oscillator constant v is varied in the 
range of 0.4 ~ 0.8 fm~ 2 . Besides the bare CD-Bonn 
and AV18 nucleon-nucleon potentials, we show the re- 
sults with using low-momentum equivalent interactions. 
The independence of the results on the oscillator con- 
stant v is demonstrated both for the ground-state energy 
and the radius. When we ignore the s\ transformation 
amplitude, the energy shows the parabolic dependence 
on v in the calculated range and the radius decreases 
monotonically with increasing v. It is interesting to see 
that the difference of the energies with and without Si is 
minimized at around the point where the radii with and 
without si coincide. 

It has been sometimes advocated |3| that different bare 
interactions collapse to universal low-momentum interac- 
tion. This property is well expected as far as the diagonal 
matrix elements are concerned, because on-shell matrix 
elements are directly related to experimental scattering 
data. However, the difference of inherent characters of 
each bare interaction which are typically characterized 
by the strength of tensor components persists in the low- 
momentum interaction. 

There are several ways to extend the present stage of 
our calculations in the future. The estimation of the 
contributions beyond the s% + s\2 approximation is nec- 
essary. The inclusion of S123, however, affects the en- 
ergy only through the change of s% and s%2, when only 
two-body interactions are present. Therefore, the net 
effect of those higher correlations may be small. The 
estimation of the contribution of a three-nucleon interac- 
tion is important to quantitatively describe the satura- 
tion properties of nuclei. The expression for the contribu- 
tion of the three-nucleon force to the ground-state energy 
is presented in Appendix D, which helps to understand 
the different character of the contributions of the three- 
nucleon force and three-body correlations (or an induced 
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three- nuclcon force). As for the treatment of short-range 
singularity of the nucleon-nucleon interaction, it is de- 
sirable to renormalize the high-momentum component of 
the bare nucleon-nucleon interaction in a more seamless 
way, working always in a harmonic oscillator basis. Fi- 
nally, numerical calculations should be extended to larger 
nuclei and also the consideration of excited states. 
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Appendix A: Similarity transformation of the one-body operator 



Remembering that the product of more than two transformation amplitudes in which the same suffix appears, e.g. 
S1S12, vanishes, the similarity-transformation of the one-body operator '^2 i U leads to up to 4-body operators. The 
expansion of the similarity-transformation e ^\ t^e with S = Sj +X)j<i s ij ~^^2i<j<k s ijk provides the following 
terms: 

1- body part ^(1 - 8^(1 + s^, (Al) 

i 

2- body part + tj)sij - Sij{U + tj) - (s. i t l + Sjtj)sij - s tJ (t t Si + tjSj)}, (A2) 

i<j 

3- body part - ^ s tj tiS ik ~ Y ^ SltiS Tjk + *.liJ' s '<- ( A3 ) 

i,j,k(j^k) i,U< k ) 

4- body part - £ (s^Usf, + s tj t iS ~), (A4) 

i,3,(.k<i) 

where the widehat ijk in the suffix means that i, j and k are to be arranged in an ascending order. 



Appendix B: Similarity transformation of the two-body interaction 

The expansion of the similarity-transformation e~ s ^X)j<j v ij^ eS provides up to 6-body operators. Explicit ex- 
pressions are the following, omitting terms including Sijki and higher excitation operators. Note that the two-body 
part corresponds to J2i<j(^ij ~ m Eq. (fTTj) . the three-body part to J2i<j<k(vijk — Uijk) originating from the 
two-body interaction Vij , and so on. 

2- body part ^(1 - Si - Sj + SiSj - sy)uy(l + S* + Sj + SjSj + s y ), (Bl) 

i<j 

3- body part ^ ~ *i ~ s j + s i s j ~ s«HjO^ + s~ + s~ S:j + s~s; + ) 

i<j k 

~( S Vk + S Tk ~ S ik Sj ~ S jk S * + S r^ V ^ ( 1 + S * + S J + S * S J + S U )}> ( B2 ) 

4- body part ^ ^{(1 - Si - Sj - s l3 + s % s 3 )v t3 (s— + + s^s— + SjS^) 

i<j k<i 

~^ S M + S W ~ SlS M ~ a i 8 $a) v ii( 1 + S * + S l + S V + s i s o)} 

+ ~ St ~ s i ~ S, J + s i s i) v ij s tk s ji + s rk s Te Vij ^ + Sl + Sj + Slj + SlSj ^ 

~( s Vk + s fk + s ^k ~ SiS Tk ~ a 3 a a) v H ( s u + s ?e + s Qt + s * s ?e + s J s ?e^> ( B3 ) 

5- body part £ ^ £{(1 - Si - Sj - Sij + s iSj )v i:j (s— + 

i<j k<l m 

+ (Srr^S-r— + S-mSr- )Vij(l + S{ + Sj + S tj + SjS,) 
v ikl jm jkl %m> Jx J J J ' 

— (s~ + s^~~ + s.-r- — SiS<~~ — SjS.~ )vij(s-rr-. + s-r~~ + SiS-rr-. + SjS<-~.) 

v im jm 13m jm J irn' J K tkt jkl jk£ J iki' 
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v jk£ jk£ J %k£ J J x vm jm ijm jm J »m /J 

/ ^ / ^ * ife jfc 0& 1 jk 3 ik' % 3 it jm 
i<Cj k^t^m 

+s u s ^ v ^ s Tk + s f k + s ^k ~ SlS Tk ~ S ' S £E )} ' (B4) 

6-body part ^ ~ s * ~ s i ~ s *i + W'*^*^ + S ikl S jm^, v v i 1 + s * + s 3 + s u + s ^i) 

z<,7 k<£ m<n 

-(s-r~~ + S-r- - SiS-r~~ - SjS^r-)Vij(s- — - + S- — - + SiS~r~~ + SjS-. — -)} 
ikK jkl jki J ikt J imn jmn jmn J imn' J 

L. J L. 1 ' J it jmn ji imn^ ik jk ijk 1 jk 3 > ik > 

i<j k.l m<n 

~{sq + S~ + S^~~ - SiS^, - SjSTljVijis^S-. + S~S- -) 

v ik jk ijk jk J ik' J v it jmn jt imn' 

+s~s~ViAs- — - + s- — - + SiS~ — - + s-jS- — -) — (s- — - + s- — - — s l s~ — - — s 1 s^—)v li s~s~ i } 

IK jt J v imn jmn jmn J imn' imn jmn jmn J imn' J ik jt J 

^ ^ ^ ^ ; ^ ik'^ jl^^^ im^ jn' ("^^) 

i<j k^l^m^n 

Appendix C: Similarity transformation of three-nucleon interaction 

A three-body interaction Vij k give the following three-body terms in the similarity-transformation. More than 
4-body operators are not shown for the sake of simplicity. 

^ (1 S{ Sj Sfc S{j &jk $ik S{Sj -(- S-iSfc ~\~ SjSfc Sijk ~1~ Si^jk Sj^ik ~t~ Sk^ij SiSjSk ~\~ S{j^ 
i<j<k 

x v ijk (l + St + Sj + s k + Sij + Sjk + Sik + SiSj + SjS fc + SjS k + s ijk + s^jk + s s ik + s k s t j + SiSjS k + s ijk ).(Cl) 

Appendix D: Expression of total energy with three-nucleon interaction 

When the Hamiltonian contains a three-nucleon interaction, X)i<j<fc v ijk, the energy obtained from Eq. (JT]) becomes 

E = J2( h M 1 + s i)\ h ) + lj2)< hh '\ Vl2 \ hh ' S > A+ 12 2x {hh'\vn\ph') A b\8i\h) 

h hh' \ p 

+ \Y, 2 x {hh'\v 12 \pp') A {p\ Sl \h){p l \ Sl \h l ) + ^J2( hh >^\PP')A(pP'\si2\hh')} 
pp> P p> 

+ 1 J2 {^^'>i23(l + si 2S )\hh'ti'} A + 3 J2( hh ' h "MP h ' h ")^\si\h) 

hh'h" p 

+ i^6 x (hh'h"\vi23\pp'h") A (p\si\h)(p'\si\h') + |j^3 x (hh'h"\vi 2 3\pp'h") A (pp'\si 2 \hh') A 



pp> pp' 



\ 9 x (htiti , \v 123 \pp'p") A (pp'\ Sl2 \hti) A (p"\ Sl \ti') \ . (Dl) 
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pp'p 



It is useful to introduce the two-body interaction ^12(3) by holding the one coordinate of 11123 with occupied states, 

(aia 2 |vi2(3) |a3a 4 )A = Y( aiCL2h \ Vl2: i\ a 3 a 4:h) A - (D2) 



h 

Using this notation, the energy Eq is written in a concise form as 



E = 5^<ft|ii(l + si)|A> + ^5^ j (hti\via + £Oi2(S)\hh'U+Y 2 x (hh'\v 12 + ^v im \ph') A (p\ Sl \h) 

h hh' I p 
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Vl2 + V 12 ( 3 )\PP') A(pp'\si2\hti) 

pp' pp' 

\ ]T \ (hh'h"\v 123 s 123 \hh'h") A + \ 9 x (htih"\v 123 \pp'p") A (pp'\s 12 \hh') A tp"\ Sl \h") \ . (D3) 



hh'h" I pp'p" 



The factors ^ and | in front of 1*12(3) in the second and third terms come from the statistical weight and are 
naturally derived. The necessity of these factors are recently pointed out by Hebeler and Schwenk in Ref. 29]. It is 
noted, however, we also have to consistently include Vijk and Sijk in decoupling equations. The expressions become 
complicated, which are beyond the scope of this article. 
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